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Abstract 

Three well known photonuclear sum rules (SR), i.e. the Thomas- Reiche-Kuhn, the 
bremsstrahlungs and the polarizability SR are calculated for 4 He with the realistic nucleon-nucleon 
potential Argonne V18 and the three-nucleon force Urbana IX. The relation between these sum 
rules and the corresponding energy weighted integrals of the cross section is discussed. Two addi- 
tional equivalences for the bremsstrahlungs SR are given, which connect it to the proton-neutron 
and neutron-neutron distances. Using them, together with our result for the bremsstrahlungs 
SR, we find a deviation from the tetrahedral symmetry of the spatial configuration of 4 He. The 
possibility to access this deviation experimentally is discussed. 
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Photonuclear sum rules (SR) are related to moments of different order of the photonuclear 
cross section and reflect important electromagnetic properties of nuclei. In fact they can 
often be expressed in terms of simple ground state properties in a model independent or 
quasi model independent way. Well known examples are the Thomas-Reiche-Kuhn (TRK) 
sum rule [3|, which gives information about the importance of exchange effects in nuclear 
dynamics via the so-called TRK enhancement factor k trk , the bremsstrahlungs sum rule 
(BSR) which is connected to the nuclear charge radius and to the mean distance 

between protons g], and the polarizability sum rule (PSR) related to the electric nuclear 
polarizability. These SR are all assuming that the dominant contribution to the cross section 
comes from unretarded electric dipole (E1UR) transitions. Two- and three-body studies 
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111 ] indeed confirm that other contributions are much smaller. Much discussed is also the 



Gerasimov-Drell-Hearn (GDH) sum rule 12], which is related to the nuclear anomalous 
magnetic moment. 

In this work we consider the TRK, BSR and PSR of 4 He within a realistic nuclear potential 
model consisting of two- and three-body forces (AV18 and UIX ^|). (The GDH sum 
rule is trivial for 4 He: it vanishes, since the 4 He total angular momentum is equal to zero.) 
We also investigate the related moments by integrating explicitly the properly weighted total 
photoabsorption cross section, which we have calculated for the same potential model jl^ . 

The aim of this study is to show that in some cases sum rules can allow to access exper- 
imentally two-body properties of the nuclear ground state, like the proton-proton, neutron- 
neutron and proton-neutron distances. In the case of 4 He this allows to test the validity of 
the configuration tetrahedral symmetry of this nucleus and at the same time "measuring" 
the amount of symmetry breaking. This work aims also at providing a guideline for experi- 
ments, where only lower bounds for the SR can be determined, as well as at giving an idea 
of the reliability of the SR approach to heavier systems, where the direct theoretical deter- 
mination of the cross section, and therefore its integration, is presently out of reach. The 
advantage to perform this kind of study in 4 He, compared to analogous ones in the two- ^| 
and three-body systems Q| , is that 4 He is a rather dense and compact nucleus, resembling 
heavier systems more closely. Only now that realistic theoretical results for the photonuclear 
cross section are available such a study is possible and one can put the extrapolation of the 
results to heavier systems on safer grounds. 

We start by recalling the formalism of the photonuclear SR. The various moments of the 
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photonuclear cross section are defined as 

m n (uj) = f duJu n af VR (u), (1) 



't»>th 

where u is the photon energy and uj t h arid u indicate threshold energy and upper integration 
limit, respectively. With indicate the unretarded dipole cross section given by 



af™{u) = QuR{u), (2) 

where Q = 4n 2 a/(3(2J + 1)) with a and Jo denoting the fine structure constant and the 
nucleus total spin, respectively The dipole response function R(lj) is given by 

R(u) = J2 IMD|0)| 2 <y(u; - E n + E ) , (3) 

n 

where \0/n) and E / n are the nuclear ground/excited state wave functions and energies, 
respectively and D is the unretarded dipole operator D = J2t=i r i T f/^i where A is the 
number of nucleons and rf and r$ are the third component of the isospin operator and the 
coordinate of the ith particle in the center of mass frame, respectively. 

Assuming that <7^ 1UR (a;) converges to zero faster than u)~ n ~ 1 and applying the closure 
property of the eigenstates of the hamiltonian H, one has the following SR for n = 0, — 1, —2: 

E TRK ee m (oo) = |(0|[D,[tf,D]]|0> (4) 
e bsr = ^(oo) = 0(o|D-D|O) (5) 
E PSR ee m_ 2 (oo) = ^^(K-^^KnlDlO)! 2 , (6) 

n 

Working out the expressions in Eqs. (JUEI), one finds that those moments are related to 
interesting properties of the system under consideration. In fact the TRK sum rule is also 
given by the well known relation 

eTRK = 6 & + kTRA ')- (7) 

where N and Z are the neutron and proton numbers, respectively, m is the nucleon mass 
and k trk is the so-called TRK enhancement factor defined as 

kTRK " |^<°I[ D ^ D P> • (8) 
From this expression it is evident that k trk embodies the exchange effects of the nuclear po- 
tential V (the double commutator in (JHJ) vanishes for systems like atoms, where no exchange 
effects are present). 



In the literature one finds a few interesting equivalences for the bremsstrahlungs sum 
rule. Rewriting the dipole operator as D = (NZ/A)Hpn, where Hpn denotes the distance 
between the proton and neutron centers of mass, one has 

S BSR =^(^) 2 (0|i?^|0). (9) 
In Ref. Foldy demonstrated that 

£ BSR =e^«>, (io) 

where (r|) is the mean square (m.s.) point proton radius 

W^(o|£rf|o>. (11) 

z i=i 

However, this relation is valid only under the assumption that the ground state wave function 
is symmetric in the space coordinates of the nucleons. 

In Ref. it was found that, in the framework of the oscillator shell model, one has 

S BSR = g ( Z 2 {r 2 } _ z<r?> ) ( (12) 

where (r^ 2 ) is the m.s. distance of protons with respect to the proton center of mass Rp 

(^}^<0|X>-Rp) 2 |0). (13) 

z i=i 

Later, in Ref. H , it was shown that the validity of Eq. (fT2|) is not limited to the oscillator 
shell model, but it is a model independent relation, which can also be written as 



^ BSR -^^ 2 W-^-^(rJ p )), (14) 
where (r 2 p ) is the m.s. proton-proton distance 

(^)^^ny(°l |>- r ;) 2 |o>- (15) 

For the BSR two additional relations exist, which are easy to prove, but which, to our 
knowledge, have not been considered in the literature, i.e. 

Z BSR = g(N*(rl)-^^(rl n )) (16) 



and 

^=g^f((rl n )-(rl)-(rl)), (17) 

where (r^) is the m.s. point neutron radius and (r£g) are the m.s. nucleon-nucleon (NN) 
distances, i.e. 

1 N 

(rL) = N(N _ 1) ( Q \ E (^i-r^lO), (18) 



N(N-1) 



Z AT 

«n> = ^(0|EE(^- r ,) 2 |0). (19) 

i=l j=l 



It is interesting to note that Eqs. (|14|) . (|16|) and (J 17)) express £ BSR via a one-body ((r„)) as 
well as a two-body quantity ((r^)). 

Finally, regarding the polarizability sum rule, one has 



Z BSR = 2n 2 a D , (20) 



where «d denotes the nuclear polarizability in the E1UR approximation. 

In this work we calculate moments and sum rules in different ways. On the one hand we 
obtain t 
section 

UIX three-nucleon force. The results have been obtained by the Lorentz integral transform 
(LIT) method Q|. The necessary equations have been solved via hyperspherical harmonics 
expansions (EIHH approach O n the other hand we obtain the SR in a more direct 



re moments by integrating our recent result for the 4 He total photoabsorption cross 
2]. This is an ab initio calculation where we use the AV18 NN potential and the 



way, as explained in the following. 

The LIT, an integral transform with a Lorentzian kernel, is defined as follows 

One way to evaluate the LIT is by using the Lanczos technique as described in Ref . |2l( . In 
fact the LIT can be re-expressed as 

£(e, r) = i (0|D • D|0) Im{(0 o |-±_|0 o )} (22) 



with z = E + e + iT and 

°|0) 



(0|D-D|0) 



(23) 
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It is evident that the LIT depends on the matrix element 

x m {z) = (0 O |__|0 O ) , (24) 

n 

which can be expressed |22( as a continued fraction containing the Lanczos coefficients 

a, = (&|#|&), h =|| b^) || , (25) 

where the <pi form the Lanczos orthonormal basis {\4>i),i = 0, . . . ,n}. Therefore the im- 
plementation of the Lanczos algorithm leads to £(e, T) (for details see jsij] ) . While the 



inversion of the LIT 23| gives access to R{uj), and thus to the moments of Eq. (0), the 
normalization of the Lanczos "pivot" \LP) = D|0) and the Lanczos coefficients allow to 
obtain the SR of Eqs. (jUinj). In fact one has: 

S PSR = Gx o(E ), (26) 
S BSR = g (LP\LP) , (27) 
S TRK = (oo - E Q ) £ bsr . (28) 

We use a HH basis, therefore the ground state |0), the Lanczos "pivot" \LP) and the Lanczos 
coefficients a n are given in terms of HH expansions. While for the ground state the expansion 
is characterized by an even hyperspherical grand-angular quantum number K and total 
isospin T=0, \LP) has to be expanded on {K' = K + 1, T = 1} states (we neglect the AV18 
isospin mixing, which is very small, as shown for the 4 He ground state in [24]). The rate 
of convergence of the various SR results from Eqs. (|26H2*B|) is given in Table |l] as a function 
of the hyperspherical grand-angular quantum number K. One observes sufficiently good 
convergence patterns for S TRK and E BSR . For the latter we have performed an additional 
test of the convergence, calculating £ BSR directly as mean value of the operator D • D on the 
ground state (Eq. (jSJ). In this way the expansion of \LP) on \K' = K + 1, T = 1} states is 
avoided. We obtain practically identical results. 

From Table |l] one sees that the convergence of S PSR is slower when the three-nucleon 



force (3NF) is included. In 15] a related problem has been found for the cross section itself. 



In fact it has been shown that the peak of the giant dipole resonance is slightly shifting 
towards lower energies with increasing K. The sum rule E PSR , which has the st rong est 
inverse energy weighting, is more sensitive to this shift than the other two SR. In [LSj an 
extrapolated cross section cr^° was obtained from a Pade approximation. We have used this 
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to determine from Eq. (JIJ the various moments for ui = 300 and 135 MeV. These results are 
also listed in Table HI 

For E TRK one sees that the SR is not yet exhausted at 300 MeV. In fact more than 20% 
of the strength is still missing. At pion threshold, one has only about 2/3 of £ TRK . As 



was already discussed in 



25| for the triton case, the rather strong contribution from higher 



energies seems to be connected to the strong short range repulsion of the AV18 potential. As 
to the TRK enhancement factor we obtain k trk = 1.31 for AV18 and 1.44 for AV18+UIX. 
These numbers are somewhat larger than older results obtained either with a variational 
wave function and AV14+UVII potential (k trk = 1.29 |26j ]) or with more approximated 
wave functions and various soft and hard core NN potentials (/t TRK = 0.9 — 1.30 (27, 28, 

For E BSR , as expected, the contribution at high energy is much smaller than for T, TRK , 
in fact at u — 300 MeV we find a missing sum rule strength of less than 2% only. For 
S PSR the strength beyond 300 MeV is even more negligible. Actually in this case the 
explicit integration leads to an even higher result than the sum rule evaluation of Eq. (|2"?j) . 
The seeming contradiction is explained by the already discussed fact that S PSR is not yet 
convergent, while for the explicit integration an extrapolated cross section is used. Indeed, 
integrating the K=18 cross section we obtain a value of 6.46 mb, which is consistent with the 
corresponding sum rule result of 6.47 mb, within the numerical error of the calculation. For 
the AV18+UIX force the value of the polarizability that we deduce from the extrapolated 
vjPSR j g Q Q555 f m s_ ^he AV18 result, which already shows a good convergence for K = 16, 
is 0.0768 fm 3 . This means that the 3NF reduces the polarizability by 15%. It would be very 
interesting to measure this nuclear polarizability by Compton scattering, as a test of the 
importance of the three-body force on such a classical low-energy observable. We find that 
vjpsr j g gj^ faofc j s affected most by the 3NF. In fact E BSR is reduced by only 10% and 
one has an opposite effect on S TRK with a 5% increase. The quenching or enhancement of 
SR due to the 3NF is the reflection of its effects on the cross section i.e a decrease of the 
peak and an increase of the tail. 

Here we would like to add a few words about the very old question, already discussed 
in ^], of the "existence" of the SR (finiteness of m n (uj) for uj — > 00), which is connected to 
the high-energy fall-off of the E1UR cross section. Since we find a rather good consistency 
between the SR and the moment values, we can state with a rather high degree of confidence 
that S TRK (and consequently the other two SR) "exist". Therefore we can try to extract 
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some information about the high energy behavior of the cross section and hence about the 
"existence" of SR with higher n. With an u~ p ansatz for the fall-off of the cross section 
above pion threshold, and requiring that £ TRK is 146 mb MeV (see Table EJ), one gets a 
rather weak energy fall-off, i.e. p ~ 1.5. This value is also consistent with £ BSR . In fact, 
adding such a tail contribution to m_i(135) one gets £ BSR = 2.39 mb, to be compared with 
2.41 mb in Table |l] The value of p might be somewhat different for other potentials, but 
probably it will not change much. Therefore we can rather safely conclude that higher order 
SR do not exist for realistic nuclear potential models. 

Finally we return to discuss the BSR, because it presents a very interesting aspect. As 
already mentioned £ BSR contains information about one- and two-body densities via (r 2 ) 
and (r^g), respectively. This means that a measurement of £ BSR and the knowledge of 
the experimental m.s. radius allow to determine (r^p) via Eqs. ([14 )1 .()16 )) and (fTTj) . In this 
way one gets information about the internal configuration of 4 He as it is explained in the 
following. 

In his derivation of Eq. (fTTl|) . Foldy assumed a totally symmetric 4 He spatial wave function, 
which corresponds to a configuration where the four nucleons are located at the four vertexes 
of a tetrahedron. For such a configuration one has (r 2 ) = (r 2 ) = (r 2 ) and (r^) = (r 2 n ) = 
( r n P ) with Qt = (r^p) I (r 2 ) = 8/3. Foldy's assumption is a very good approximation for 4 He, 
but other spatial symmetries (mixed symmetry, antisymmetric) are also possible. What can 
be learned from our S B result with respect to this question? For 4 He, which is a T=0 
system one can safely assume that (r 2 ) = (r 2 ) = (r 2 ) (isospin mi xing is tiny Using in 

Eqs. (UH), (HEJ) and (JT2I) the AV18+UIX value of (r 2 ) = 2.04 fm 2 |lj (which coincides with 
the experimental one |30j], corrected for the proton charge radius) and £ BSR from Table U we 
obtain (r 2 p ) = (r 2 n ) = 5.67 fm 2 and (r 2 n ) = 5.34 fm 2 , i.e. two values which differ by about 
6%. The ratios Q pp ( nn ) = ( r p P („„))/( r2 ) an d Qn P = (rl p }/(r 2 } are not much different from 
the correspondent value Qt of a classical tetrahedral configuration. We obtain Q pp = 2.78 
and Q np = 2.62 instead of Qt = 2.67. One notices that Q pp ( nn ) — Qt — 2(<5„, p — Qt)- This 
reflects the different numbers of proton-proton and neutron-neutron pairs (2) with respect to 
proton-neutron pairs (4). Using Eq. Q one can also derive the distance between the proton 
and neutron centers of mass. One has R pn = 1.58 fm instead of 1.65 fm for the tetrahedral 
configuration. 

Notice that with (r 2 ) = 2.04 fm 2 the "tetrahedral" BSR would be 2.62 mb, the same value 
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that one obtains using Eq. (|10|). This value is 9% larger than our result. The distortions 
that we find from our 9% smaller BSR are the consequence of the different effects of the 
potential on isospin triplet and isospin singlet pairs. 

We conclude that when considered in its body frame 4 He should look like a slightly 
deformed tetrahedron. Of course this statement has to be intended in a quantum mechanical 
sense, regarding the mean square values of the nucleon-nucleon distances on the two-body 
density. It is clear that one cannot measure this deformation "directly", since it is not a 
deformation of the one-body charge density (the 4 He charge density has only a monopole). 
On the other hand such a deformation is accessible experimentally in an indirect way via 
the measurements of the charge radius and of X BSR . 

This leads to the question how exactly £ BSR can be measured in a photonuclear experi- 
ment. Two points have to be addressed: i) the contributions of El retardation and higher 
multipoles, which are not contained in E BSR , but which will contribute to the experimental 
cross section and ii) the contribution of the high energy tail. To this end it is instructive to 
consider the results from the three- nucleon photodisintegration. Regarding point i), addi- 
tional effects of the El retardation and higher multipoles have been calculated in Ref. 3| 
for the AV18+UIX potentials in a Faddeev calculation. Using those results one finds that 
El retardation and higher multipoles increase m_i(135) by about 1% only. In fact according 
to Gerasimov[^| there is a large cancellation of El retardation and other multipoles. There 
is no reason for a larger effect in 4 He. On the contrary, since the leading isovector magnetic 
dipole (Ml) transition is suppressed in this nucleus (it is zero for an S-wave) one can expect 
an even smaller contribution. As to point ii) considering the fall-off of the triton E1UR 

n 

cross section around pion threshold from Ref. |32(] one obtains a result very similar to our 
4 He case, namely p ~ 1.5. However, including the other multipole contributions one gets a 
considerably smaller value, namely p ~ 1.1. This is no contradiction with the small increase 
of 1% for E B , since the inverse energy weighted cross section integrated from 100 to 135 
MeV gives only a rather small contribution to the sum rule. On the other hand one would 
overestimate £ BSR taking the full cross section. Thus we suggest that the tail contribution 
to an experimental BSR should be estimated using the theoretically established fall-off oj~ p 
with p ~ 1.5. 

In conclusion, we have studied three well known photonuclear sum rules for 4 He within 
a realistic two- and three-nucleon potential model. Two new equivalences for the BSR sum 



9 



rules have allowed us to deduce information about two-body properties of the nuclear ground 
state, like the proton-proton, neutron-neutron and proton-neutron distances. In particular 
we have tested the validity of the configuration tetrahedral symmetry of this nucleus and 
have found that this symmetry is slightly broken. We have suggested an experimental way 
to access this symmetry breaking via a measurement of the BSR which could be performed 
in one of the existing or planned low-intermediate energy photonuclear facility. 
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TABLE I: Convergence in K of the SR for AV18+UIX potentials. The converged AV18 results 
are also shown. The last two lines of the table show the convergence in Q of the various moments. 





AV18+UIX 




S PSR 


VjBSR 


yjTRK 


K 


10~ 2 [mb MeV- 1 ] 


[mb] 


10 2 [mb MeV] 


8 


6.230 


2.398 


1.430 


10 


6.277 


2.396 


1.448 


12 


6.331 


2.394 


1.451 


14 


6.382 


2.401 


1.458 


16 


6.434 


2.406 


1.460 


1 Q 


6.473 


2.410 


1.462 




AV18 




7.681 


2.696 


1.383 




AV18+UIX 


uj [MeV] 




m_i (lo) 


m (co) 


135 


6.55 


2.27 


.944 


300 


6.55 


2.37 


1.14 
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